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Determination of low energy parameters for NN-scattering at N 4 LO in all partial 

waves with j ' < 5. 
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The Variable S-matrix approach offers a unique way to extract low energy threshold parameters 
for a given NN potential. We extract those parameters for the np system from the Nijmll and 
Reid93 potentials, to all partial waves with total angular momentum j < 5, using the generalized 
effective range expansion 
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where f SJ = (S SJ — l) /(2ifc) is the scattering amplitude and S SJ is the unitary S-matrix in coupled 
channel space with total spin s and total angular momentum j. Our calculation includes all the 
relevant contributions of the full amplitude to order 0(k 8 ) in the CM momentum. We also discuss 
the validity of the generalized effective range expansion in the region of analyticity k < m w /2. 

PACS numbers: 03.65.Nk,13.75.Cs,21.30.Fe 
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I. INTRODUCTION 



The study of the NN interaction has played a dominant 
role in the theory of nuclear structure y|. At low ener- 
gies the corresponding scattering phase shifts can be best 
parameterized in terms of an effective range expansion 
(ERE) 0, 0| ■ This expansion can be suitably generalized 
for all partial waves and coupled channels (for a review 
see e.g. Ref. Q and references therein.). The calculation 
of the low energy threshold parameters is straightforward 
in potential scattering, and can be extracted from the 
asymptotic form of the wave function computed order by 
order in a low energy expansion, but unlike potentials 
the ERE parameters are not subjected to off-shell ambi- 
guities. Benchmark descriptions of NN scattering data 
have been obtained by the phase shift partial wave anal- 
ysis (PWA) of the Nijmegen group carried out a decade 
ago p| and further parameterizations through high qual- 
ity potentials |(|. For channels involving central waves 
low energy threshold parameters have been determined 
routinely and for the high quality potentials this calcula- 
tion has been undertaken in Ref. [?J). However, a detailed 
and systematic determination of these threshold param- 
eters in non central partial waves is, to our knowledge, 
still missing. 

Almost parallel with the previous developments there 
has been in the last decade a renewed interest on the 
NN interaction based on the application of Lagrangian 
effective field theories (EFT) methods, motivated by 
Weinberg's ideas Q | B and pursued by many oth- 
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El El 111 111 El US EH!! SIM El to attempt a real- 
istic description of np scattering data. This approach has 
open promising perspectives in the theory of nuclear sys- 
tems (for a recent review see e.g. Ref. 136| and references 
therein). Actually, the reference database most modern 
EFT calculations are confronted with are those of the Ni- 
jmegen group [10]. Moreover, the determination of the 
low energy constants (LEC's) appearing in the EFT La- 
grangian depends on the regularization scale and also on 
the short distance behaviour of the long distance chiral 
potential due to One, Two and higher Pion exchanges. 
The results for the scattering matrix and in particular for 
the low energy threshold parameters should, of course, be 
independent on the renormalization prescription. In the 
absence of explicit pion exchanges the EFT appr oach can 
be mapped into an effective range expansion |13| , and the 
LEC's can, in principle, be deduced from the physical 
threshold parameters for a given renormalization scale. 
When long distance chiral potentials generated by suc- 
cessive pion exchanges are included the situation becomes 
much more involved, particularly if those potentials can- 
not be treated perturbatively. Recently, we have pro- 
posed a regularization scheme |34l El where the LEC's 
can also be deduced directly from the low energy thresh- 
old parameters beyond perturbation theory within a vari- 
able S-matrix approach Is?! . The implementation of such 
an approach requires an explicit knowledge of the thresh- 
old parameters as input. In practice, most EFT calcu- 
lations fit their results to the Nijmegen database H H 
in a given energy window in terms of the LEC's within 
a given regularization scheme and at a given scale, and 
low energy threshold parameters could be determined af- 
terwards. Under these circumstances we think it of in- 
terest to provide in any case these low energy threshold 
parameters as deduced directly from the high quality po- 
tentials HH in all partial waves. 
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In this paper we determine the low energy parameters 
of the Nijmll and Reid93 potentials 5] for all partial 
waves with j < 5. To this end we use the variable S- 
matrix ap pro ach (for a review with many references see 
e.g. Ref. |37]]. An application of the variable phase ap- 
proach in the context of EFT in the singlet 1 Sq channel 
can be found in Ref. |22j) where the S-matrix can be 
computed by solving a non-linear system of first order 
differential equations representing the variation of the S- 
matrix under a continuous switching on of the potential 
starting from a trivial potential (see below for a sketchy 
derivation). By implementing a low energy expansion 
one obtains a set of first order coupled differential equa- 
tions from which the low energy threshold parameters 
can be obtained as asymptotic values at long distances 
for any partial wave. The method is stable and whenever 
comparison can be made reproduces previous known re- 
sults for the ERE parameters in central waves. Partial 
aspects of the present paper have already appeared 
in connection with the 1 S'o and 3 Si — 3 D\ channels and 
in the context of the renormalization of the singular OPE 
potential. 

The paper is organized as follows. In Sect. [H] we in- 
troduce our notation for scattering amplitudes, partial 
wave expansion M-matrix and scaled M-matrix, which 
enable to take the low energy expansion in a particularly 
transparent and clean way. We also discuss the interplay 
between low energy expansion, rotational covariance and 
unitarity. In Sect. IIIII we briefly outline the derivation 
of the variable S-matrix for coupled channels, following 
the method of adiabatically switching on the potential 
from the origin. Translating this equation to the variable 
scaled M-matrix we write down in Sect. ITvl the coupled 
differential equations which determine as a final condition 
the low energy threshold parameters. Our numerical re- 
sults are presented in Sect.[V] In Sect. IV Al we show the 
deduced threshold parameters for the Nijmll and Reid93 
potentials for all partial waves with total angular momen- 
tum j < 5. In Sect. IV 51 we examine the validity of the 
expansion for these partial waves. In Sect. IV CI we study 
the poles and residues of the scattering matrix, mainly in 
connection with the 1 5'o virtual state and the 3 Si — 3 D\ 
bound state deuteron within the coupled channel effec- 
tive range expansion. Finally, in Sect. lVll we present some 
conclusions and outlook for further work. In Appendix lAl 
we provide some raison d'etre for the present calculation 
and show the difficulties we have encountered when a di- 
rect fit to the database is attempted. In the Appendix FBI 
we investigate further the virtual and bound state prop- 
erties in an expansion of the inverse scattering length in 
the s— wave channel. 



that for the two nucleon system in a total spin s and 
total isospin t state, the elastic scattering amplitude for 
the transition k, m s — > k , m' s with fixed CM energy E = 
k 2 /2[i n p reads 1 

f^, ms (fc',fc) = X) C 'fi'(*'.*)[ 1 -(-) ,+i+t ]^ (!) 
jlV 

where 

Cfj,(fc',fc) = yV~' C(l,s,j\m - m a ,m s ,m)Yi tfn - ma (k) 

m 

x C(l',s,j\m-m' s ,m' s ,m)Yi>, m - m , a (k')* (2) 

and f ; S j, is the scattering amplitude for the 2s+1 lj state 
with a total angular momentum j, total spin s and orbital 
angular momentum I and is given by 



fSj _ _}_ ( QSj 



5u> 



(3) 



in terms of the symmetric and unitary S-matrix, Sf, 



St 3 ,,, Y<v S*\iS$ v — 5u>. Due to parity conservation, 
for j > the states with I = j cannot couple to states 
with I = j ±1. Thus, for the spin singlet state, s = 0, 
one has I — j and hence the state is uncoupled 



(4) 



whereas for the spin triplet state s = 1, one has the 
uncoupled I — j state 



S 



1.7 



e j 



(5) 



and the two channel coupled states 1,1' — j ± I states 
which we use Stapp-Ypsilantis-Metropolis (SYM or Nu- 
clear bar) |4d| | parameterization for 



qij _ ( - s / : ( I S j-1]+1 | 

( fc,- \ 2iS lj , • • (c- \ i(S l3 ,+S 1 i,) 

_ I cos(2ej)e j- 1 ism(2eyje n i- 1 
\i sin (2ej)e y j- 1 cos(2ej)e 3+ 1 



In the discussion of low energy properties it is also inter- 
esting to use the Blatt-Biedenharn (BB or Eigen phase) 
parameterization |4l) that is given by 



gij _ 1 cosfj — smej 



sine,- 



cos e, 



COS Co- 



sine.; 



2ia lj ' , 
e ^ 



a. 

e j +1 



(6) 



II. S-MATRIX, PHASE SHIFTS AND 

THRESHOLD PARAMETERS 1 ^ e use ^ e normalization for the differential cross section 



As it is well known (see e.g. Ref. [3|f ), the requirement 
of all relevant symmetries in the NN-interaction implies 



d*(k,m s ^k',m> s ) = t . - 
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The relation between the BB and SYM phase shifts is 
given by 



5 lj +S lj - S lj +5 lj 



tan(2e,-) 
tan(2ej) 



(7) 
(8) 



In order to take the low energy limit we define the 
standard symmetric and real coupled channel scaled M- 
matrix 



S = (M + ifcD 2 ) (M - ifcD 2 ) 



(9) 



with k the CM np momentum and D = 
diag(fc Zl , . . . ,k N ). Due to unitarity of the S-matrix in 
the low energy limit, k — ► we have 



(S - 1)„ , = -2ia v ,ik 



i'+i+i 



(10) 



with ayi the (hermitian) scattering length matrix. The 
threshold behaviour acquires its simplest form in the 
SYM representation, 



- 




(11) 




- -aV^ +1 , 


(12) 


% - 




(13) 






(14) 




4 -tt«fc 2 i+i. 


(15) 



defining a natural hierarchy for phase-shifts (see also be- 
low) 2 . The scaled M- matrix, M, has a good low energy 
behaviour Q| and admits the coupled channel analog of 
the effective range expansion 



M = a 



-rk 2 

2 



v 2 fc 4 + v 3 fc 6 + v 4 A: 8 



, (16) 



where a, r and v 2 are the scattering length matrix, 
effective range and curvature parameters respectively. 
Higher order matrix parameters will be denoted for mu- 
tational simplicity V3, V4, and so on. Within the ERE 
our notation is as follows: LO means including order k° 
terms, NLO including k 2 terms and NNLO including fc 4 
terms, and so on. The expansion in k 2 around the ori- 
gin holds up to the next singularity, which in the case 
of NN interaction corresponds to the pion left cut, so 
\k\ < ±777,7,. /2 ~ 70MeV. One of the advantages of the 
low energy expansion at the level of the scaled M-matrix 
is that unitarity is preserved exactly at any order in the 
expansion. 

If one diagonalizes the M-matrix using the orthogonal 
transformation, which we call fi, in Eq. © 



m = DOD M DO D 



(17) 



2 In the BB form one has similar behaviours for the 5's but Ej 



the corresponding eigenvalues are related to the BB eigen 
phaseshifts as follows 



m; = k 2l+1 cotd° J 



. 1 +1 fk 2 +v?k 4 + ... ,(18) 
a, £ 



for I = j ± 1. 

In connection with the low energy expansion note that 
the total amplitude Q satisfies both the off-forward op- 
tical theorem from the unitarity of the S-matrix (dq rep- 
resents solid angle interaction in the unit vector q direc- 
tion), 

^m' g m s ( ) ) ^m s m' s V s ) = 
dqC Jm X k '^)CjrnX k ^r (19) 

m'J 

as well as rotational covariance, i.e. under rotations k — » 
Rk one has 3 

m'J ,m'J' 

x C sms (Rk\Rk)D s m ,,, trns (R)* 

(20) 

with D s m , m //(R) the rotation Wigner matrices corre- 
sponding to spin s. This induces a mixing of the cou- 
pled channel states, since they are irreducible under ro- 
tations. The latter holds for any truncated sum in the 
partial wave expansion in Eq. (Q. On the other hand, 
the low energy expansion of the amplitude only satisfies 
this properties perturbatively; the coefficients of k and 
k' in a Taylor expansion of the amplitude f^, m (k,k) 
are neither rotational covariance nor do they satisfy ex- 
act unitarity. This can be seen from the coupled channel 
amplitude, 



f = D M - ikD z 



D 



(21) 



which yields after taking the k — > expansion yields 

1 



DaD -DaraDfc 2 



(22) 



This threshold behaviour sets up a natural low energy 
hierarchy for the 2s+1 lj states 

O(k ) 1 S a 3 S 1 

0{k 2 ) 1 P 1 3 P Q 3 P 1 3 P 2 ,E 1 

0(k 4 ) 1 D 1 3 D 1 3 D 2 3 D 3 ,E 2 

0(k 6 ) 1 F 3 3 F 2 3 F 3 3 F 4 ,E 3 



0(k 8 ) 1 G 4 3 G 3 3 G 4 3 G 5 ,E 4 



(23) 



Corresponding to the invariance of f(k',k). 
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where we use the notation E\ 



Si - 3 Di , E 2 



P 2 ~ 3 F 2 , E 3 = 3 D 3 - 3 G 3 , E A = 3 F 4 - 3 H A , and so on. 
(in general Ej = 3 (j + l)j — 3 (j — l)j). This means that 
for a complete calculation at 0(k 8 ) we have to include 
the a's for the states 1 G*4 , 3 G3 , 3 G4 , 3 G5 and^, a's and 
r's for 1 F 3 , 3 F 2 , 3 F 3 , 3 and E 3 , a's, r's and V2's for 
1 £>i , 3 £>i , 3 £>2 , 3 D3 and £2 and a's, r's V2's and i>3's for 
1 Pi , 3 P , 3 -Pi , 3 P2,Ei, and a's, r's, i^'s, v 3 's and v^s for 
1 5'o and 3 Si. In general, to 0(k 2n ) one needs 



AT 



2 x (n + 1) + 5 x -n(n + 1) 



(24) 



independent parameters. Thus, at O(k ) we need only 
N = 2 parameters, at C(fc 2 ) one has N = 9 parameters, 
at C(fc 4 ) , JV = 21, at C(/c 6 ) , JV = 38 and at C(fc 8 ) 
we have JV = 60 parameters, etc. Although these seem 
rather large numbers, note that the declared number of 
independent adjustable parameters to NN data in the 
high quality potential of Ref. Q is 41. The very recent 
EFT calculation NNNLO in the Weinberg counting H3 
declares 26 adjustable parameters (actually LEC's) to the 
NN database of Ref. . This calculation corresponds to 
go to order 0(k 6 ) for the theory without explicit pions, 
i.e. for k <C m^/2. 

In this paper we use the low energy expansion in a way 
that both unitarity (|19fl and rotational covariance 1)20)) 
are simultaneously satisfied. As a consequence, the low 
energy expansion for the amplitude is not complete in 
the sense that starting at a given order on k not all 
the contributions of higher order are taken into account. 
Nevertheless, up to order C(fc 8 ) we compute all the 60 
parameters. The highest angular momentum involving 
this order is j = 5. For completeness we present all the 
effective range parameters up to V4 in all partial waves 
with j < 5. We will see below, however, that presenting 
the data in terms of the scaled M-matrix provides some 
interesting insights regarding the validity of the effective 
range expansion. 



III. 



VARIABLE S-MATRIX 



In this Section we re-derive the variable S-matrix equa- 
tion using a continuous deformation of the potential. Al- 
ternative derivations are based on Jost functions |37j and 
inner boundary conditions |.3 ll l35| . The coupled channel 
Schrodinger equation for the relative motion reads 

r l 2 ' 

U(r) + 3 



u"(r) 



u(r) = k u(r) , 



(25) 



where U(r) = 2/x n?) V(r) is the coupled channel matrix 
potential which can be written as for j > 0, 



(26) 



U lj (r 





C#(r) 




u 



We will take for these potentials the ones available in 
Ref. 0. InEq. \ 2 = diag(h(h + 1), . . . ,l N (l N + 1)) 
is the angular momentum, u(r) is the reduced matrix 
wave function and k the CM. momentum. In the case at 
hand N = 1 for the spin singlet channel with / = j and 
JV = 3 for the spin triplet channel with l\ = j — 1, I2 = j 
and I3 = j + 1. For ease of notation we will keep the 
compact matrix notation of Eq. 1)25(1. The potentials in 
Ref. are regular at the origin, so the regular solution 
is given by the boundary condition at the origin 



u(0) = 0, 



(27) 



At long distances, we assume the asymptotic normaliza- 
tion condition 



u(r) -> h ( - } (r) -h (+) (r)S 



(28) 



with S the standard coupled channel unitary S-matrix. 
The corresponding out-going and in-going free spherical 
waves are given by 



h (±) W = di ag (hf i (kr),...,hl(kr)), 



(29) 



with hif(x) the reduced Hankel functions of order I, 
hf{x) — xH^^^x) ( h = e ±lx ), and satisfy the free 
Schrodinger's equation for a free particle. 

In order to deduce a variable .S-matrix equation, we 
determine first the infinitesimal change of the S matrix 
under a general deformation of the potential U(r) — > 
U(r)+AU(r). Using Schrodinger's equation l|25|) and the 
standard Lagrange's identity adapted to this particular 
case, we get 

[u(r)tAu'(r) - u'(r)t Au(r)] ' = u(r)'AU(r)u( r ) 

(30) 

which, after integration from the origin to infinity and 
using the asymptotic form of the matrix wave function, 
Eq. 1)28)1 . as well as the regular condition at the origin, 
Eq. (23 yields 



2i/«S t AS=/ dru(r)^AU(r)u(r 
Jo 



(31) 



In particular, for the parametric family of potentials 
U(r, R) = 6{R - r)U(r) we get AU(r, R) = S(R - 
r)U(r)Ai? and hence 



2ifcS t (i?)S'(i?) = u(R)^V(R)n(R) , 



(32) 



and using the value of the wave function at the outer 
boundary, Eq. ffijjl 



u(R) = h^(R)-h^(R)S{R), 
we finally get the variable S-matrix equation, 



(33) 



2ifc 



dS(R) 
dR 



S(R)h i+) {R) -&-\R)] U{R) 



h^(R)-h^ + \R)S(R) 



(34) 
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This is a first order non-linear matrix differential equa- 
tion which can be solved by standard means, provided the 
S-matrix is known at one given scale. Note that for any 
value of the boundary radius we have a different on-shell 
scattering problem. In the case of a regular potential, 
Eq. (|34|) has to be supplemented with an initial condi- 
tion at the origin, namely the trivial one (corresponding 
to the absence of a potential), and its asymptotic value 
yields the full S— matrix ; 



S(0) = 1 



S = S(oo) 



(35) 



IV. EVOLUTION OF LOW ENERGY 
PARAMETERS 

In order to take the low energy limit of Eq. 1)34(1 and 
corrections there-off, we introduce the variable or running 
M— matrix, in analogy with Eq. Q 



S(R) = (M(R) + ifcD 2 J (M(R) - ifcD 2 J , (36) 
as well as the reduced Bessel functions 

j = g (h« - h<")) , (37) 

-y = i(h<+> + h<->) , (38) 

with 

jl[x) = xji(x) , yi(x) = xyi(x) . (39) 



Thus, we get 

M'(k,R) = ^M(i?,fc)ij(fci?)D" 1 -y(kR)nj U(R) 

x (pikR^^MiR, k) - y(kR)Tij . (40) 

The scaled variable M-matrix admits the analog of the 
effective range expansion 

M(R) = -a(R)- 1 + ^r(R)k 2 + v 2 (i?)fc 4 + v 3 (i?)fc 6 + . . .(,41) 

where a(i?), r(i?), v 2 (i?), v 3 (_R), v 4 (i?) , .. are the cor- 
responding running parameters. In this form the low 
energy limit can be easily taken. Defining the matrix 
functions 

B k (R) = diag (y h (kR)k h ,...,yi N {kR)k l ») (42) 
and their low energy expansion 



A k (R) = J-i^D- 1 = A + fc 2 A 2 + fc 4 A 4 + . . . . 
k 

B k (R) = y(fci?)D - B + fc 2 B 2 + fc 4 B 4 + . . . , 



(43) 



we get the system of coupled equations 



— [a^R)]- 1 = -([a(i?)]- 1 A + B )U(i?)(Ao[a( J R)]- 1 +B ) , 

■^r(R) = ([a(i?)]- 1 A + B )U(i?)(r(i?)A + 2[a(i?)]- 1 A 2 + 2B 2 ) 

+ (r(R)A + 2[a(R)]- 1 A 2 + 2B 2 ) U(R) ([a(R)}- l A + B ) , (44) 

^v 2 (i?) = ([a(i?)]- 1 A +B )U(i?) ^-[a(i?)]- 1 A 4 + ir(i?)A 2 + v 2 (i?)A -B 4 j 

+ ^-[a(J?)]- 1 A 4 + \v(R)A 2 + v 2 (i?)A - B 4 ^j U(R) ([a(i?)]- 1 A + B ) 

+ ( lr(R)A 2 - [a(R)]- 1 A 2 - B 2 ) U(R) (l-r(R)A 2 - [a(i?)]- x A 2 - B 2 ) , 



and similar equations for V3(i?) and v 4 (i?). These equa- 
tions generalize to the coupled channel case those already 
found in Ref. [3^| and have to be supplemented with the 
initial conditions, 



The physical threshold parameters correspond to the val- 
ues at infinite, 



a(0)=0, r(0) = 0, v 2 (0) = (45) 



a(oo)=a, r(oo) = r, v 2 (oo)=v 2 ,...(46) 
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In the triplet coupled channel case the SYM threshold 
parameters matrices are 



where we get for the effective range parameters the rela- 









tions, 


a J-l 


*n 


(47) 


«3S1 
^3S1 


r j 


*n 
#J' 


(48) 


«3S1 


v i-l 

V 




(49) 





and so on. Similar definitions hold for the uncoupled 
channels. 



V. NUMERICAL RESULTS 
A. Determination of low energy parameters 

Low energy scattering data can directly be described 
in terms of threshold parameters, like a, r, etc. , de- 
fined through Eq. l(Tr]|l . Unfortunately, besides a and ro 
in the 1 Sq and 3 Si channels, the partial wave analysis 
data base p| does not provide values for them except 
for V2,Va,V4 in Ref. in the deuteron channel. In Ap- 
pendix^] we show that a direct fit to the database turns 
out to be numerically unreliable. The NN data base pro- 
vides explicit potentials like the Nijmll and Reid93 po- 
tentials, for which the variable phase approach may di- 
rectly be applied. In such a way we can uniquely and 
accurately determine all the needed low energy threshold 
parameters by integrating Eqs. Ij44(l from the origin to 
infinity with the trivial initial conditions Eq. (|45|) . In all 
cases we have checked that the variable S-matrix equa- 
tion, Eq. (JH1J, indeed reproduces the results of the NN 
data base |5j for the phase shifts. In practice, our results 
are stable if one integrates up to Roo = 40fm. Our results 
for the low energy threshold parameters for the Nijmll 
and Reid93 potentials can be summarized in Table [I] To 
obtain these results we looked for a stable plateau in the 
large R region in each threshold parameter separately. So 
we quote in Tabic [I] the stable digits in this flat region. 
This is the main result of this work. In general we see 
a rather good agreement between both potentials never 
worse than 5% and frequently much better. We also ob- 
serve, as expected, that numerical accuracy worsens by 
going to higher orders in the ERE in a given partial wave 
or going to higher partial waves. 

In the interesting case of the 3 Si eigen- channel, one 
has 4 



kcotS 3S i 



^— + Irrssik 2 + vssik 4 + . . . (50) 

C*3S1 * 



2r E \OiE\ r 3D1 a 2 E1 



«3S1 



(51) 



X 3S1 



1 



A - E1 (2a 3D1 r 3D1 r E1 - a E1 r E1 + 8v m ) 
4a 3 Si 

_ 2 

. -f 1 (a 3 Dirl D1 - 2a E if 3 DifEi + 4u 3 m) 

4a 3Sl 



1 



4 «3S1 



(52) 



and so on. Using the numerical values for the Nijmll 
(Rcid93) potentials from Table |U we get 

a 3S i = 5.41896 (5.42293) fm 

r 3S1 = 1.75334 (1.75556) fm 

v 3S i = 0.04531 (0.03266) fm 3 

v' 3S1 = 0.65831 (0.65797) fm 5 

V ^ S1 = -4.19144 (-4.19262) fm 7 



(53) 



in agreement with previous findings 0, El • 



4 We use the notation v = V2,v' = v% and v" = V4 for simplicity 



B. Validity of the effective range expansion for 
partial waves 

Once the low energy threshold parameters have been 
determined we can readily check the validity of the ef- 
fective range expansion given by Eq. (|16|) . On theoreti- 
cal grounds the full expansion around the origin should 
be convergent within the analyticity domain, i.e. up 
to the next singularity, which for NN interaction hap- 
pens at about k = i.e. Slab ~ — lOMeV 
where it develops a logarithmic branch cut due to the 
One Pion Exchange Potential. At k = iim^ a Two 
Pion Exchange dilogarithmic branch cut sets in at about 
-Slab ~ — 40MeV, and so on. From this viewpoint any 
realistic description of the phase shifts should be consid- 
ered unreliable for energies beyond the analyticity do- 
main ~ lOMeV. On the other hand, this does not nec- 
essarily imply that the polylogarithmic corrections are 
numerically large. So one may directly learn from the 
"data" how important these corrections are in practical 
terms. There are two possible ways to do this, either 
by using the SYM phase shifts given in Ref. 5] or by 
directly evaluation of the scaled M-matrix defined. In 
terms of the LAB energy, Slab = k 2 /fi np the scaled 
M matrix can be represented as a constant, straight line 
and a parabola if we keep the scattering lengths, effective 
ranges and curvature parameters. The results of such a 
comparison for the scaled M-matrix is given in the series 
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TABLE I: Low energy threshold parameters for all partial waves in np scattering for the Nijmll and Reid93 (in brackets) 
potentials. We give the scattering length a, the effective range r and the curvature parameters V2, Vs and V4,. Units are in 
relevant powers of a fm, with I ,1 = j,j ± 1. From them one can compute the M- matrix at low energies using M i ; /fc k — 
-(a.- 1 ) l}V + i(r) M /fc 2 + (v 2 ) M 'fc 4 + (v 3 ) u ,fc 6 + (v 4 ) M /fc 8 + . . . where M = (S - f ) (S + ik and S is the unitary S-matrix in 
coupled channel space. In the coupled channel case we use the SYM(Nuclear bar) low energy parameters. All presented digits 
are numerically significant. A long dash ( — ) stands for cases where numerical accuracy was outraged and no value could be 
reliably deduced. 



Wave 


a(fm i+i '+ 1 ) 


ro(fm !+i ' +1 ) 


w 2 (fm i+r + 3 ) 




^(fm i+r+7 ) 




Nijmll ( Reid93 ) 


Nijmll ( Reid93 ) 


Nijmll ( Reid93 ) 


Nijmll ( Reid93 ) 


Nijmll ( Reid93 ) 


L S 


-23.727(-23.735) 


2.670(2.753) 


-0.4759(-0.4942) 


3.962(3.652) 


-19.88(-18.30) 


3 Po 


-2.468(-2.469) 


3.914(3.870) 


1.099(0.9616) 


3.816(3.712) 


-7.6(-7.4) 


l Pi 


2.797(2.736) 


-6.399(-6.606) 


-1.580(-1.834) 


0.404(1.024) 


7.8(8.4) 


3 Pi 


1.529(1.530) 


-8.580(-8.556) 


-0.0181(+0.01006) 


-0.86(-0.934) 


0.3(0.2) 


"S 1 


5.418(5.422) 


1.833(1.833) 


-0.131(-0.141) 


1.444(1.433) 


-7.95(-7.9) 


3 D 1 


6.505(6.453) 


-3.523(-3.566) 


-3.699(-3.803) 


1.12(1.023) 


-3(-2.6) 


Ei 


1.647(1.645) 


0.404(0.413) 


-0.274(-0.264) 


1.447(1.423) 


-7.3(-7.3) 


L D 2 


-1.389(-1.377) 


14.87(15.04) 


16.37(16.73) 


-13.2(-13.0) 


34.0(34.0) 


*D 2 


-7.405(-7.411) 


2.858(2.851) 


2.395(2.370) 


-0.99(-1.00) 


2.0(2.0) 


3 P 2 


-0.2844(-0.2892) 


-8.270(-8.363) 


-6.91(-7.13) 


-6.0(-6.3) 


-20(-30) 


3 F 2 


-0.9763(-0.9698) 


-5.640(-5.821) 


-22.96(-23.78) 


-79.32(-83.0) 


-117(-127) 


E 2 


1.609(1.600) 


-15.70(-15.89) 


-25.18(-25.72) 


-23.3(-24.8) 


-67(-70) 


l F 3 


8.383(8.365) 


-3.924(-3.936) 


-9.888(-9.937) 


-15.4(-15.6) 


-4(-4) 




2.703(2.686) 


-9.932(-9.994) 


-20.56 (-20.73) 


-19(-19.4) 


-20(-20) 




-0.1449(-0.1770) 


1.369(1.365) 


2.06(2.040) 


2(2) 


-(-) 


3 G 3 


4.880(4.874) 


-0.03306(-0.0529) 


-0.0166(-0.1261) 


-0.117(-0.662) 


-3.1(-4.1) 


E 3 


-9.695(-9.683) 


3.255(3.249) 


7.655(7.618) 


9.5(9.3) 


-(-) 


L G, 


-3.229(-3.210) 


10.78(10.82) 


34.4(34.50) 


80(84) 


-(-) 


3 G 4 


-19.17(-19.14) 


2.056(2.059) 


6.810 (6.822) 


16.6(16.7) 


10(14) 


3 F 4 


-0.01045(-0.01053) 


-3.02(-3.04) 


-8.0(-8.0) 


-10(-) 


-(-) 


3 m 


-1.250(-1.240) 


-0.1762(-0.178) 


-1.52(-1.527) 


-ll(-ll) 


-35 (-40) 


Ei 


3.609(3.586) 


-9.489(-9.539) 


-29.60(-29.72) 


-68(-69) 


-80(-70) 


L H 5 


28.61(28.57) 


-1.724(-1.727) 


-7.925(-7.92) 


-32(-32.5) 


-60(-60) 




6.128(6.082) 


-6.41(-6.45) 


-24.9(-25.1) 


-87(-86) 


-(-) 


y G 5 


-0.0090(-0.010) 


0.48(0.481) 


1.9 (1.874) 


6(6) 


-(-) 


3 h 


10.68(10.66) 


0.0108(0.0107) 


0.145(0.143) 10" 5 


1.43(1.4) 


6.4(6) 


E 5 


-31.34(-31.29) 


1.554(1.553) 


6.99 (7.018) 


28(28) 


50(50) 



of figures Fig. JTJ), @, 10, I0J, © and ® for all partial 
waves with < j < 5. In table |n] we present the effective 
range parameters for each partial wave in units of powers 
of m n /2, so that they correspond to the contribution of 
different terms of the ERE at the maximum value of k 
within the domain of analyticity of M(fc). 

On the light of these figures, two conclusions may be 
drawn. In the first place, for a fixed LAB energy as 
we increase the angular momentum, the relative error to 
the full scaled M-matrix increases. This is expected since 
higher angular momenta are more sensitive to the long 
distance physics, which is pion dominated. In the limit 
of large angular momentum peripheral waves, the low 
energy threshold parameters should be described solely 
in terms of pion dynamics. The second observation is 
that the absolute error decreases for increasing angular 
momentum as it corresponds for a suppression of large 
scattering angles. 



C. Poles and Residues of the scattering matrix 



The scattering amplitude has poles for negative ener- 
gies 5 



7 + ik 



(54) 



where 7 > corresponds to a bound state (first Riemann 
sheet in E) and 7 < (second Riemann sheet in E) to 
a virtual state. The coefficients Af J correspond to the 
asymptotic bound state wave function 



(55) 



We assume that there are no degenerate poles. In the two channel 
case this condition is automatically satisfied. 
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TABLE II: Contributions Mo, Mi, Ma, M3 and M 4 for all partial waves in np scattering for the Nijmll and Reid93 (in 
brackets) potentials at a center of mass momentum of ' m n /2 (equivalent to a laboratory energy about 10 MeV). From them 
one can compute the M-matrix at fc cm = m w /2 using M;.;/ = (Mo);.;' + (Mi)^/ + (Ma);.;' + (M3);.;/ + (M 4 );.;/ + . . . where 
M = D (S — 1) (S + 1) _1 ifcD, S is the unitary S-matrix in coupled channel space and D = diag(fc ix , . . . , fc ijv ) is the coupled 
channel centrifugal factor. M = -a -1 , Mi = ±r (m^/2) 2 , M 2 = v 2 (m^/2) 4 , M 3 = v 3 (m^/2) 6 , M 4 = v 4 (m^/2) 8 , and so 
on. We have taken m w — 138.0MeV. A long dash ( — ) stands for cases where numerical accuracy was outraged and no value 
could be reliably deduced. 



Wave 


M 


Mi 


M 2 


M 3 


M 4 




Nijmll ( Reid93 ) 


Nijmll ( Reid93 ) 


Nijmll ( Reid93 ) 


Nijmll ( Reid93 ) 


Nijmll ( Reid93 ) 


L S 


0.0422(0.0421) 


0.1633(0.1683) 


-0.0071(-0.0074) 


0.0072(0.0067) 


-0.0044(-0.0041) 


3 Po 


0.4052(0.4050) 


0.2353(0.2366) 


0.0164(0.0144) 


0.0070(0.0068) 


-0.0017(-0.0017) 


± 1 


-0 3575C-0 36551 


-0 3912(-0 40391 


-0 0236(-0 02741 


0007(0 00191 


0017(0 00191 


3 Pi 


-u.ujoy^-u.uooo j 




-u.uuuo^-u.uuuz J 


n nni nm 7^ 

-u.uuiu^-u.uui 1 J 


U.UUUl ^U.UUUUO ) 


3 5i 


-0.1999(-0.1999) 


0.1121(0.1121) 


-0.0020(-0.0021) 


0.0027(0.0026) 


-0.0017(-0.0017) 




-0.1666 (-0.1680) 


-0.2154(-0.2180) 


-0.0553(-0.0569) 


0.0021(0.0019) 


-0.0006(-0.0006) 


£1 


0.05062(0.05097) 


0.02467(0.02524) 


-0.00410(-0.00395) 


0.00266(0.00261) 


-0.00164(-0.0016) 


is A 


720(0 7261 


909(0 91 91 


245(0 2501 


-0 024(-0 0241 


008(0 0081 


3 Do 


1 350(0 1 3491 


1 747(0 1 7431 


0358(0 03541 


-0 001 8(-0 001 81 


0005(0 00051 


3 P 2 


-0.4222(-0.4257) 


-0.5056(-0.511) 


-0.1033(-0.1065) 


-0.011(-0.012) 


-0.005(-0.01) 


3 F 2 


-0.1230(-0.127) 


-0.3448(-0.356) 


-0.3433(-0.356) 


-0.1450(-0.152) 


-0.026(-0.029) 


E 2 


-0.6960(-0.7022) 


-0.9601 (-0.9712) 


-0.3764(-0.3845) 


-0.043(-0.045) 


-0.015(-0.02) 




-0.1193(-0.1196) 


-0.2399(-0.2406) 


-0.1478(-0.1486) 


-0.0281(-0.0285) 


-0.001(-0.001) 


*F* 


u.ouyy^ u.o; ^oj 


n 6079( n 6i 1 nl 




o n^fi( n n^6l 




^ 3 


0.05153(0.05151) 


0.0837(0.08347) 


0.0308(0.0304) 


0.004(0.004) 


-(-) 




-1.530(-1.870) 10~ 3 


-2.021(-3.233) 10~ 3 


-0.248(- 1.885) 10" 3 


-0.21(-1.21) 10" 3 


-0.7(-0.9) 10" 3 


E 3 


0.1024(0.1023) 


0.1990(0.1986) 


0.1144(0.1139) 


0.0174(0.0170) 


-(-) 


L G 4 


0.3087(0.3115) 


0.6587(0.6616) 


0.514(0.516) 


0.15(0.15) 


0.02(— ) 


3 G 4 


0.0522(0.0522) 


0.1257(0.1259) 


0.1018(0.1020) 


0.0304(0.0305) 


0.003(0.003) 


3 F 4 


-0.0961(-0.0965) 


-0.185(-0.186) 


-0.12(-0.12) 


-0.02(— ) 


-(-) 


3 H 4 


-0.00080(-0.00082) 


-0.01077(-0.01088) 


-0.0227(-0.0228) 


-0.020(-0.020) 


-0.01(-0.007) 


Ei 


-0.2774(-0.2791) 


-0.5801 (-0.5831) 


-0.443(-0.445) 


-0.125(-0.125) 


-0.02(-0.02) 


L H 5 


-0.0350(-0.0350) 


-0.1054(-0.1055) 


-0.118(-0.119) 


-0.059(-0.059) 


-0.015(-0.02) 




-0.1632 (-0.1644) 


-0.392(-0.394) 


-0.373(-0.38) 


-0.16(-0.15) 


-(-) 


3 G 5 


0.0109(0.0109) 


0.0294(0.0294) 


0.028(0.028) 


0.01(0.012) 


-(-) 


3 T 

J-5 


-0.0092(-0.0105) 10" 3 


0.661(0.65) 10" 3 


2.2(2.2) 10" 3 


2.6(2.6) 10" 3 


1.5(1.4) 10 -3 


E 5 


0.0319(0.0319) 


0.0950(0.0951) 


0.1045(0.1046) 


0.051(0.051) 


0.01(0.011) 



whereas the relativistic expression is given by 

Bf = M p + M n - yjM* - ^ - sjMZ-rt (61) 

The calculation of threshold properties in the 3 5i — 3 D\ 
channel from given deuteron properties is a standard pro- 
cedure (see e.g. Ref. 0). Here we do just the opposite, 
i.e. compute deuteron properties from threshold parame- 
ters within the the effective range expansion. In practice 
this means solving the equation 

Det(M(fc)- 1 - ikT> 2 )\ k= , n = (62) 

using the LO,NLO and NNLO approximations. Residues 
are evaluated by numerical integration using Cauchy's 
theorem. The numerical results for the smallest poles in 
the singlet 1 5 , and triplet 3 5i — 3 Di channels and their 
corresponding residues are presented in Table ITlTl As can 



In the 1 So channel one has 7 = —j v and the asymptotic 
wave function is 

uiso(r) -> Ai S0 e lvr (56) 

In the case of the deuteron one has 7 = 7d 

«ssi(r) -> A 3S1 e~^ r (57) 

3 

The ratio s to d wave is defined as 



u 3D i(r) ^ A 3D1 e-™ ( 1 + — + ) (58) 



A3D1 /_„s 

The non-relativistic deuteron binding energy reads 

B? R = £- (60) 
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be seen the expansion is convergent but even to N 4 LO in 
the amplitude is not sufficient to reproduce the observ- 
ables in a completely satisfactory way within experimen- 
tal uncertainties. The difference should in principle be 
attributed to higher orders in the momentum expansion 
as well as other effects. In appendix |B] we analyze these 
results further on the light of an expansion around the 
limit of large s-wave scattering lengths. Our conclusion 
is that such an expansion does not improve on the de- 
scription of the bound state observables. 



VI. CONCLUSIONS AND OUTLOOK 

In this paper we have extracted the low energy thresh- 
old parameters for all np partial waves up to states with 
total angular momentum j < 5 taking into account the 
coupled channel nature of the problem for realistic NN 
potentials. Our description entails up to order 0(k 8 ) 
in the CM momentum of the full np scattering ampli- 
tude. An adequate framework is to use the effective 
range expansion of the scaled coupled channel M-matrix, 
where the kinematic centrifugal factors have been fac- 
torized out. These low energy parameters comprise the 
scattering length matrix a, the effective range matrix r, 
the curvature matrix V2 and higher like V3 and V4 and 
are relevant for a shape independent description of the 
scattering data in terms of the scaled M-matrix in a re- 
gion of analyticity in the complex energy plane around 
the origin which radius extends up to the left partial 
wave logarithmic cut generated by One Pion Exchange 
intermediate states. The practical determination of these 
parameters from the solutions of the Schrodinger equa- 
tion may be cumbersome, so we have found extremely 
convenient to use the variable S-matrix formalism. The 
low energy threshold parameters can be directly deter- 
mined from the asymptotic solution of a set of coupled 
non linear differential equations which correspond to an 
adiabatic switching on of the NN potential using trivial 
initial conditions. Finally, we have found that the cou- 
pled channel effective range expansion works well within 
the expected region of analyticity, namely k < and 
in fact a clear trend to convergence is observed. However, 
beyond the region of analyticity we do not expect the low 
energy expansion to be realistic regardless on how many 
terms are included in the expansion. Actually, the ana- 
lyticity domain of the amplitude can be enlarged if the 
proper left cut singularities, corresponding to OPE,TPE, 
etc. are implemented. For such a program the method 
presented in our previous work |34l l35l | looks particularly 
promising. Finally, one problem with the use of high 
quality potentials has to do with the determination of 
errors on the potentials and hence on the low energy pa- 
rameters. In this paper we have used two such potentials 
to assess those errors, but it would be rather interesting 
to make the error analysis directly based on the coupled 
channel effective range expansion. 
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APPENDIX A: DIFFICULTIES IN EXTRACTING 
THE LOW ENERGY PARAMETERS FROM A 
FIT 

At first sight one might think that the effective range 
parameters could be determined directly from a fit to 
the Nijmegen data base ,5] and avoid the use of the 
corresponding potentials. In this appendix we want to 
elaborate on the problems we have encountered while fit- 
ting that data base within a generalized coupled chan- 
nel effective range expansion of Eq. . Unfortunately, 
this data base does not provide error estimates for their 
phase shifts (although 8 significant digits are given), nor 
the typical energy resolution where these data should be 
trusted, so some compromise must be made. 

We use the NN-data and define the % 2 as 

where we take A£lab = O.OIMeV , and M NN and AM NN 
are the mean value and the standard deviation of the six 
potentials listed in the NN-data base [5j, which can be 
taken as independent uncorrelated primary data. The 
factor M/(4k) is the Jacobian of the transformation be- 
tween the Lab-energy and the CM. momentum, £lab = 
2k 2 /M, and would correspond to make an equidistant 
sampling in p, in the limit AE'lab — > (this is why we 
take a small energy spacing) . This weight factor is intro- 
duced in order to enhance the region at low momenta. 
On the other hand, very low momenta must be excluded 
since the resulting mean value M-matrix is incompatible 
within the attributed errors with the expected theoret- 
ical behaviour, Eq. (| 1 €>l) . so we take -Elab > 0.5MeV. 
This is partly due to the poor accuracy of the data at 
low energies; the calculation of the scaled M-matrix re- 
quires increasing accuracy at low energies. Also, the fit 
goes up to £lab < lOMeV, which corresponds to a CM. 
momentum about k — m n /2 where we expect the finite 
polynomial of the scaled M — matrix to truly represent 
an analytical function within the convergence radius up 
to the branch cut singularity located at k = ±im OT /2. 
The form of the fitting function is 

Mer = v q + vik 2 + v 2 k 4 + v 3 k 6 + v 4 k s + ... (A2) 

In Fig. [7] we show as an illustration the 1*2 parameter 
determined from a fit to the low energy region of the 
NN data base as a function of the maximal LAB- 
energy considered in the fit. As we see, instead of a 
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TABLE III: Poles and residues of the scattering amplitude for the Nijmll and Reid93 (in brackets ) p otentials based on the 
effective range expansion, Eq. I|16|l . The data for Bd are from |42j, for rjd from |43j and for As from |44|. 





LO 

Nijmll ( Reid93 ) 


NLO 
Nijmll ( Reid93 ) 


NNLO 
Nijmll ( Reid93 ) 


N 3 LO 
Nijmll ( Reid93 ) 


N 4 LO 
Nijmll ( Reid93 ) 


Exp. 




-0.04215(-0.04213) 


-0.04001(-0.03994) 


-0.04001 (-0.03994) 


-0.04001(-0.03994) 


-0.04001 (-0.03994) 




B% R (MeV) 
Bf(MeV) 

A&nifnT 1 ' 2 ) 
Ai D1 {frn-V 2 ) 


0.18472(0.18438) 
1.41197(1.40990) 
1.41250(1.41043) 
0.01036(0.01033) 
0.60732(0.60701) 
0.00629(0.00627) 


0.22511(0.22498) 
2.10146(2.09909) 
2.10263(2.10026) 
0.02669(0.02659) 
0.79351(0.79330) 
0.02117(0.02110) 


0.23050(0.23045) 
2.20345(2.20251) 
2.20475(2.20380) 
0.02512(0.02502) 
0.86620(0.86711) 
0.02176(0.02170) 


0.23127(0.23120) 
2.21807(2.21684) 
2.21938(2.21815) 
0.02519(0.02510) 
0.87805(0.87873) 
0.02212(0.02205) 


0.23148(0.23141) 
2.22215(2.22088) 
2.22347(2.22219) 
0.02521(0.02512) 
0.88219(0.88283) 
0.02224(0.02218) 


2.224575(9) 
2.224575(9) 
0.0256(4) 
0.8846(9) 



plateau within some energy window, we observe an ever 
changing value. We observe no stability depending on the 
number of terms considered in Eq. (|A2|) either. For com- 
parison we also plot the values we obtained by integrat- 
ing the Eqs. fHJl with the Nijmll potential in Sect. ITVl 
which where quite stable numerically. As we see, the val- 
ues obtained from the fit, in the chosen energy window 
are hardly compatible. The deceptive features extend to 
other channels, and non diagonal low energy threshold 
parameters such as the matrix elements of a and r. 

Finally, we have also tried, with no success, other 
methods for the determination of the low energy thresh- 
old parameters, like evaluation of derivatives within sev- 
eral algorithms. The reason for the failure has to do with 
round-off errors generated by the relatively small number 
of digits provided in the NN database. Actually, at very 
low energies these round-off errors make the construction 
of the scaled M matrix itself rather unstable numerically, 
since large centrifugal factors l/k l are involved. 



APPENDIX B: THE LIMIT OF LARGE 
SCATTERING LENGTHS 

In this appendix we discuss some interesting aspects of 
the limit of large scattering lengths. For a s— wave eigen 
phase shift, i.e. 1 S'o and 3 Si eigen channels, the effective 
range expansion can be written as 

fccot<5 = -- + \r Q k 2 + v 2 k i + v 3 k 6 + v 4 k 8 + . . .(Bl) 
ao 2 

The poles in the S-matrix are given by the solutions of 
the equation 



/o (£7) 1 = k cot S — ik 







k— vy 



(B2) 



For a given order of the truncated effective range ex- 
pansion the previous equation reduces to an algebraic 
equation with real coefficients and it has at least one real 
solution, which can be taken as the analytical continu- 
ation of the lower order ones by taking the limit of the 
effective range parameters rp, v 2 , ■ ■ • — > 0, or alternatively 



as the limit ao — > 00. All other roots either wander to 
infinity or accumulate to build, in the limit of infinitely 
many terms of the expansion, a branch cut. Expanding 
into powers of I/ckq one gets 



1 



7 



^o 
2al 



5 r 3 



V 2 



O' 



2al 



«3 



( 



33 6 
16 r 



15rgv 2 



4 V 2 



4r i>3) 



+ 



(B3) 



The neglected terms contain the parameters v±, . . .. 
Within the same approximation the residue for the pole 
of the S-matrix in this (eigen)channel defined in analogy 
to Eq. (53J, Sq -> Al/(~f d + ik) becomes 



, 2 2 3r hr 2 
a at ag 



M r g - 42r v 2 2|ir5 - 126^ + 2v 3 



- 330v 2 r% + 12v 3 r + 12v\ 
a 7 



(B4) 



The expansions Eq. I|B3|) and Eq. (|B4|) are convergent 
provided ao is within the domain of analyticity of the 
exact solution of Eq. I|B2|) . Otherwise the expansion is 
asymptotic, i.e. we can use it to evaluate numerically 
the value of 7 in the limit ao — * 00, up to a given order 
where the remainder starts increasing basically due to the 
presence of large factorials. Nevertheless one can use the 
expansion with a given error estimate. For instance, at 
NLO the domain of analyticity is given by the condition 
of a vanishing discriminant of a second order algebraic 
equation yielding ao = 2tq = 3.50652 for the analyticity 
inner boundary. The numerical values for ao and tq lie 
within the boundary and one can use the series expansion 
to any order to evaluate 7 with increasing accuracy. In 
the NNLO case the analyticity boundary is given by an 
analogous discriminant condition. Numerically we find 
the inner boundaries for ao located at lower points ao = 
3.444, -0.777, -0.133 for the Nijmll potential. 
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TABLE IV: Poles and residues of the scattering amplitude for the Nijmll potential based on the effective range expansion, in 
the large scattering length limit a — » oo to C(1/q 5 ) with a — ai Sg in the 1 So channel, a = «3si in the 3 Si — 3 Di channel ( 3 5i 
eigen channel), and a = 0351 in the "bar" 3 Si — 3 Di channel (barred quantities). This calculation is complete to N 4 LO (see 
main text). 





Oil/a) 
Nijmll ( Reid93 ) 


0(l/a 2 ) 
Nijmll ( Reid93 ) 


0(l/o?) 
Nijmll ( Reid93 ) 


0(l/a 4 ) 
Nijmll ( Reid93 ) 


0(l/a 5 ) 
Nijmll ( Reid93 ) 


T^fm- 1 ) 


-0.04215(-0.4213) 


-0.03977(-0.03969) 


-0.04004(-0.03997) 


-0.4000(-0.03997) 


-0.4001(-0.03994) 


7d(fm~ i ) 
B? R (MeV) 
B R (MeV) 


0.18454(0.18440) 

0.0(0.0) 

0.0(0.0) 

0.0(0.0) 
0.36908(0.36880) 


0.21439(0.21425) 
1.41227(1.41020) 
1.41227(1.41020) 

0.0(0.0) 
0.54820(0.54789) 


0.22405(0.22391) 
1.86922(1.86673) 
1.86922(1.86673) 
0.01035(0.01032) 
0.64479(0.64452) 


0.22791(0.22779) 
2.05403(2.05146) 
2.05456(2.05199) 
0.01218(0.01212) 
0.69896(0.69889) 


0.22962(0.22952) 
2.13694(2.13461) 
2.13782(2.13549) 
0.01671(0.01663) 
0.73010(0.73027) 


7"d(fm _1 ) 
B% R {MeV) 
Bf(MeV) 
f]d 


0.18454(0.18440) 

0.0(0.0) 

0.0(0.0) 

0.0(0.0) 
0.36907(0.36880) 


0.21576(0.21558) 
1.41227(1.41020) 
1.41227(1.41020) 

0.0(0.0) 
0.55639(0.55584) 


0.23050(0.23037) 
1.89010(1.88689) 
1.89010(1.88689) 
0.01035(0.01032) 
0.68707(0.68679) 


0.23209(0.23197) 
2.15612(2.15363) 
2.15665(2.15416) 
0.01233(0.01227) 
0.74308(0.74290) 


0.23413(0.23404) 
2.21874(2.21629) 
2.21963(2.21718) 
0.01741(0.01733) 
0.78694(0.78718) 



TABLE V: Same as table|^but for higher orders. 





0(l/a e ) 
Nijmll ( Reid93 ) 


0(l/a 7 ) 
Nijmll ( Reid93 ) 


0(l/a 8 ) 
Nijmll ( Reid93 ) 


0(l/a 9 ) 
Nijmll ( Reid93 ) 


N 4 LO 
Nijmll ( Reid93 ) 


T^fm" 1 ) 


-0.04001 (-0.03994) 


-0.04001(-0.03994) 


-0.04001(-0.03994) 


-0.04001(-0.03994) 


-0.04001(-0.3994) 


7d(fnT 1 ) 
B% R (MeV) 
Bf(MeV) 
Vd 


0.23047(0.23038) 
2.17666(2.17462) 
2.17773(2.17568) 
0.01848(0.01838) 
0.74867(0.74907) 


0.23091(0.23083) 
2.19696(2.19515) 
2.19814(2.19632) 
0.02053(0.02042) 
0.76000(0.76059) 


0.23116(0.23109) 
2.20784(2.20620) 
2.20908(2.20744) 
0.02167(0.02156) 
0.76717(0.76790) 


0.23131(0.23123) 
2.21396(2.21245) 
2.21523(2.21371) 
0.01953(0.01942) 
0.77181(0.77265) 


0.23148(0.23141) 
2.22215(2.22088) 
2.22347(2.22219) 
0.02521(0.02512) 
0.77826(0.77939) 


7"d(fm _1 ) 
B? R (MeV) 
B R (UeV) 
Vd 

(Atsrfifm- 1 ) 


0.23344(0.23336) 
2.26308(2.26123) 
2.26423(2.26238) 
0.01895(0.01886) 
0.79951(0.79996) 


0.23356(0.23349) 
2.25980(2.25807) 
2.26106(2.25934) 
0.02149(0.02140) 
0.81118(0.81198) 


0.23289(0.23282) 
2.26231(2.26093) 
2.26365(2.26226) 
0.02238(0.02228) 
0.80941(0.81033) 


0.23272(0.23266) 
2.25179(2.25044) 
2.25314(2.25179) 
0.01979(0.01969) 
0.80990(0.81102) 


0.23148(0.23141) 
2.22215(2.22088) 
2.22347(2.22219) 
0.02521(0.02512) 
0.77826(0.77939) 



If we instead make the large scattering length expan- 
sion in the determinant bound state condition, Eq. (|62ll . 
for a truncated scaled M-matrix we get an algebraic equa- 
tion in the deuteron pole k = vy which order is deter- 
mined by the level of approximation. In the SYM pa- 
rameterization the LO, NLO, NNLO, NNNLO effective 
range parameters start contributing at order 7, 7 2 , 7 4 
and 7 6 respectively and also contribute to a maximal or- 
der 7 6 , 7 7 , 7 9 and 7 12 respectively. This means that 
the coefficients of 7 in an expansion in 1/0:351 are com- 
plete to 0(03"^) at LO, O^^) at NLO, and CKa^) 
at NNLO. This can also be seen by applying the large 
scattering length expansion of Ea. lB3l to the eigen phase 
shift, Eq. (|50|> . with the identifications i|52|) . yielding for 
the pole 

1 . 1 r 3S i . 1 1-2 . - - \ — 3 

1 "351 1 

1 / 2 5 , n- - 

+ 7 ( ^OL E1 r iD i + -r 3S1 + 8a E ir 3 sirEi 



-a. 3D ifEi-^v 3 ,sija~g 1 + ... (B5) 
We also get for the residues 

A d A d 25l ei gsoirgi - 5a E iflsi , m fi s 

A 3S1 A 3D1 = — T2 1" ■ • • (Bfaj 

"3S1 a 3Sl 

Finally, for the mixing coefficient rjd we get 

Vd = — + ^4 + •■■ (B7) 

a 3Sl za 3Sl 

and so on. Note that the leading contributions in the in- 
verse scattering length behave differently for any observ- 
able. According to Eq. fB"^ . making a 3 si large makes 
also di3si large. Thus one may study separately the con- 
vergence of observables in both cases. The numerical 
values for the pole position and corresponding residua 
evaluated in the inverse scattering length expansion are 
presented in Tables llVl for the expansion up to order 1 /a 5 
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and 1 /a 5 . This order involves for pole position only up to 
N 4 LO in the amplitude (for the binding energy and mix- 
ing parameter the situation is somewhat different, see 
e.g. Eq. 1B7I) ). As we have already discussed the con- 
verse is not true, and N 4 LO contains a series of higher- 
power corrections in the inverse scattering length. We 
present those corresponding only to N 4 LO in Table [V] 



up to ninth order in the inverse scattering length. As 
we see, and one could have anticipated, the convergence 
for the virtual state pole is faster than for the deuteron 
bound state. Actually, we see that the expansion in the 
inverse scattering length for bound state properties is less 
convergent than the ERE pursued in Sect. IV 01 
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Lab. Energy ( MeV) 



Lab. Energy ( MeV) 



FIG. 1: np scaled M-matrix on the effective range expansion for states with total angular momentum j — 0. We construct the 
Nijmll scales M-matrix corresponding to the generalized effective range expansion of Eq. Ijlfip to a given order as a function 
of the LAB energy. LO means including order k° terms, NLO including k 2 terms and NNLO including k 4 terms with the low 
energy threshold parameters obtained from solving the evolution equations for the threshold parameters, Eq. I|44^ with the 
Nijmll potential . Data are obtained from the partial waves phase shift analysis from Ref. pj- 
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FIG. 2: np scaled M-matrix based on the effective range expansion for total angular momentum j = 1. For notation see Fig.Q 
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FIG. 3: np scaled M-matrix based on the effective range expansion for total angular momentum j = 2. For notation see Fig.Q 
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FIG. 4: np scaled M-matrix based on the effective range expansion for total angular momentum j = 3. For notation see Fig.Q 
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ER-NNLO 
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FIG. 5: np scaled M-matrix based on the effective range expansion for total angular momentum j = 4. For notation see Fig.Q 
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FIG. 6: np scaled M-matrix based on the effective range expansion for total angular momentum j = 5. For notation see Fig.Q 
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5 6 7 8 9 10 11 12 13 14 15 5 6 7 8 9 10 11 12 13 14 15 

Ei^ Max (MeV) E, ab Max (MeV) 



FIG. 7: The ^parameter for the 1 So (left) and the 3 Si (right) channels determined from a fit to the low energy date of the 
NN data base [3 (see Eq. IjAl^ l and main text, as a function of the maximal LAB-energy considered in the fit. p n means a 
fit including up to 0(k n ) terms in the effective range expansion Eg . IA2I "Database" means a fit to the average value of the 
corresponding scaled M-matrix. "Reid93" and "Nijmll" means a fit to only this data. The values we obtained by integrating 
the Eqs. (glj Nijmll and Reid93 potentials. 



